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ABSTRACT 

Assuming that perturbative QCD is the dominant explanation for the narrowness of 
the vector quarkonia, we perform a x 2 minimization analysis of their hadronic decays as 
a function of two parameters, the mass of the gluino and the value of as(Mz). A value 
below 1 GeV for the gluino mass is strongly preferred. Consequences for SUSY breaking 
scenarios are discussed. 
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Recently it has been pointed out that the quarkonia data can be made consistent with 
minimal supersymmetric grand unification if the gluino and photino masses are below half 
the Z° mass [[TJ . In addition, the quarkonia data becomes consistent with the world average 
measurements of the strong fine structure constant, a^(Mz), if these masses are not above 
the T region. The result of Ref. |3J was based on differences in the running of the strong 
coupling constant in the presence of gluinos. If, however the gluino mass is below half the 
T mass a second effect comes into play Namely, the possibility for a quarkonium state to 
decay into gluino-containing final states affects the analysis of as at the relevant scale for 
that decay In the current work we seek a best fit value to the gluino mass, m g , taking both 
effects into account. This allows us to consider the full range of relevant gluino masses in 
a x 2 analysis and to determine the best fit value of the gluino mass m g . 

This work depends on the following basic assumption: The dominant explanation for 
the narrowness of the quarkonia states including the J/tj) and <fi is perturbative QCD. 

Any future models for the quarkonia data that rely on very large non-perturbative 
effects or relativistic corrections should be considered as alternative to the present analysis 
and judged on the basis of their relative physical plausibility. Small non-perturbative 
corrections, of course, would only modify our results by small amounts. 

We analyze up to second order in QCD the hadronic decay rates of six quarkonia states, 
0(1019), J/^(3097), J/^(3686), T(9460), T(10020), and T(10350). Ignoring possible 
higher order or non-perturbative effects and relativistic corrections each of these defines 
within errors the strong coupling constant at an appropriate scale 

«3 OwsC*)) = a 3 ,i ± Si, (i = l,...,6). (1) 

Assuming the gluino mass, m g , lies above half the mass of the quarkonium state, these 
values are independent of the gluino mass and are given in Ref. using the 1992 branching 
ratio averages of the Particle Data Group @. In the current work we consider also lower 
values of m g so that the az,i become functions of m g . 

If the gluino mass lies below half the quarkonium mass, there are decays into two 
gluinos (gg) and into two gluinos plus one gluon (ggg) which compete with the standard 
three gluon (ggg) decay. These gluino-containing final states || however are suppressed by 
four powers of m q /mq and are negligible for currently allowed values of the squark mass rriq. 
The dominant gluino-related correction to the quarkonium width is therefore the two-gluon 
plus two-gluino (gggg) final state with no intermediate squarks. The corresponding decay 
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rate as a function of the gluino to quark mass ratio m g /m q was written down in Ref. 
by applying a color +spin correction factor to the rate for the ggqq decay as calculated in 
Ref. §. The result is 

r( 3 Si(gg) -» gggg) = 3a 3 (fis,mg) 
T( 3 S 1 (qq) ^ ggg) n 

where r = 2m g /m( 3 Si) ~ m g /m q . For small r (but r > 0.1), 
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32(tt 2 - 9) 
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i?(r) = -InrH p-£ r - 3.56 - r 2 (8 lnr + 0.95) + — ttV 
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1 6 

+ r 4 ( ln 2 r + 1.1 lnr -8.1) + C(r 5 ) . (3) 



Eq. (0) is infrared divergent as r — > 0. A correct treatment || shows that -R(O) ~ 1.57. 
We use the above expression for R in region r > 0.1 and use a quadratic interpolating 
polynomial to join R smoothly to its value at r = 0. The results presented here are not 
sensitive to the exact form of the cutoff. R(r) falls rapidly for increasing r. For r > 0.5 it 
is necessary to use the exact integral expression for R(r). It is possible |lj to choose the 
scale ns(i) so that the known (first order) standard model corrections to the 3 gluon decay 
rate vanish identically. The hadronic decay rate into states of dissimilar quarks is then 

rfSxiqq) - hadrons) = r{ 3 Si(qq) - ggg) + TfS^qq) - gggg) + ■■■ (4) 

We have then 

3/ \(, , 3a 3 ,i(^S,^ff) D/ 3 /r\ 

a 3)l (^s, mj) (1 + y -R{r) J = a 3>i . (5) 

The gggg and gggg correction to Eq. ([|) is absorbed into the three gluon term by the 
choice of scale (is Q- The coupling on the right hand side of Eq. (|j) is the result 
obtained in Ref. [Q assuming no gluino contribution {m g > m( 3 Si)/2). 

Our analysis proceeds as follows. For fixed value of the gluino mass we solve Eq. (|5|) 
to obtain six values of the strong coupling constant at six scales appropriate to each of the 
vector quarkonia. The extrapolation to these scales from the Z° mass is done using the 
three loop renormalization group equation, 

^da 3 (Q) <A(Q) r b . _ N c 2/ ^l 
Q— B? 1 = 7^ a + 4^ a 3(Q) + 16^2 "1(3) > ( 6 ) 



dQ 2tv 



where 



° = 11 -3 n ^( 1+ loS + 2S)- 2 ^ (7a) 

6 = 102-yn / (Q)-48n § (Q), (76) 

2857 5033 325 2 ,_. 

^8"^^) + — nJ(Q). (7c) 

The effects of scalar quarks and Higgs bosons are neglected in accord with the decoupling 
theorem. The off-diagonal two loop effects are treated as electroweak corrections to the 
one-loop a coefficient. We neglect the running of the electroweak corrections, using instead 
as an average value 

a em = 1/133, sin 2 #vK = 0.2333, (8) 



and 



5 Ctem Q^em / r .\ 

3 cos^tV sm 9 W 



The gluino contribution to the 3-loop c coefficient is not known. We keep the zero gluino 
contribution (Eq. ([7c])) but monitor the effect of the c term as an estimate of unknown 
perturbative effects. In the fits of reasonable x 2 , the contribution of the c term is small. 
We take into account the threshold dependence of quarks and gluinos according to the 
formula 

Q'' 



n * (Q) = / (J&)' (106) 



with 



f{x) = 1 + . 1 In 

2Jx(l + x" 



y/l + X — ^fx 



:id 



-y/l + X + y/x- 

For each value of a^[Mz) and nig we can extrapolate from the Z° mass down to the 
quarkonium region using Eq. and calculate a x 2 ■ 

The basic assumption given above does not require the total absence of non- 
perturbative effects. In fact, in the case of the (most accurately measured) IS" quarkonia 
states the graphs shown in Ref. 0] do exhibit some scatter in the data around the QCD 
predictions at about the 10% level. To find the preferred values of the gluino mass we 
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follow two alternative customary procedures for the treatment of data in such cases. The 
results of these two procedures do not greatly differ. 

Method 1: discard the data point of worst agreement and minimize the % 2 of the 
remaining data points. In the present case this is the J/ip(lS) decay. The remaining five 
vector quarkonia decays, which contain at least one entry from strange, charm, and bottom 
quarks agree well with the theory and provide a relatively sharp minimum x 2 as a function 
of two parameters a^{Mz) and the gluino mass rrig. 

Method 2: retain all six vector quarkonia decays but add in quadrature with the 
experimental errors a "theoretical error" to take into account possible non-perturbative or 
binding effects, i.e. 

<*3(j*s(i)) = <X3,i(ji S , m~g) ± sj5f + Xjal ^s^rrig). (12) 

On general grounds one would expect such corrections to be more important for the lighter 
quarkonia than for the T states. The A^ parameterize our ignorance about higher order 
and non-perturbative corrections. They do not, of course, constitute a model for such 
effects (since they do not shift the central values of 013^) and in fact no reliable model 
exists apart from lattice QCD which has not as yet attained sufficient numerical accuracy. 
Clearly for sufficiently large Xi all predictive power is lost. We take our basic assumption 
to imply Xi <C 1. An adequate % 2 is found with = 0.05 for the bottomonium states, 
and Xi = 0.10 for the charmonium and strangeonium vector states. Our conclusions are 
qualitatively insensitive to increasing this value for the lighter quarkonia in the sense that 
the gluino mass of minimum x 2 remains low although the x 2 values increase more slowly 
away from the minimum. Similarly the favored light gluino also persists for smaller Xi 
although the minimum x 2 is then not a mathematically acceptable fit. 

For comparison with the best fits allowing a light gluino, we show in Fig. 1 the best 
fit to the quarkonia data in method 2 assuming the gluino lies at high mass (400 GeV) 
where it essentially decouples. This fit seems surprisingly good to the naked eye specially 
considering that it relates pure perturbative QCD to the hadronic decay rates of six vector 
quarkonia of three species over mass scales varying by a factor of ten. However it is not 
a mathematically good fit since it corresponds to a x 2 per degree of freedom (x 2 /DoF) of 
3.7. In Fig. 2 we show the variation in the x 2 /6 for this heavy gluino case as a function of 
as(Mz). The minimum x 2 of 3.7 is several standard deviations worse than the best fits 
with a light gluino. In addition this best fit corresponds to a value of a 3 (M^) that is many 
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standard deviations away from the world average value and is inconsistent with SUSY 
unification with a SUSY scale below lOTeV. A similar attempt to fit the five quarkonia 
states as in method 1 but without light gluinos would yield a minimum x 2 P er degree of 
freedom many times larger than the minimum x 2 of Fig. 2. 

In Fig. 3 we show the x 2 contours for method 1 treating the gluino mass rrig and 
a 3 (M z ) as variable. The x 2 /5 = 1 contours define two acceptable regions: 

a 3 (M z ) = 0.1135 ±0.0005, m- g = (0.32 ± 0.05) GeV; (13a) 



or 



a 3 (M z ) = 0-1145+^, m~ g = (o.Ol +°^) GeV. (136) 



The best fit corresponding to the central values of ( |13a ) is shown in Fig. 4. Comparing with 



Fig. 1, it is clear that the fit has improved due to decay of T states into gluino containing 
hadrons and due to the slower falloff of the coupling constant as a function of energy. In the 
solution corresponding to ( |136| ), even the <p decay has significant contribution from gluino 
containing final states. In this case both gluinos must presumably hadronize into a single 
pion where they can readily mix with gluon pairs. With light gluinos one must expect 
that all hadrons have non-negligible gluino components just as there is a non- negligible 
probability to find strange quarks in the sea of non-strange hadrons. 

Fig. 5 shows the x 2 contours for method 2. In this method the x 2 /6 = 1 contour lies 
within the region: 

a 3 (M z ) = 0.1115 ^o'ooJg, m g = (0-44 ± 0.17) GeV. (14) 

More conservative values (90% CL) can be read from the x 2 /F)oF = 2 contours in Fig. 3 
or Fig. 5. In method 2 there is also a tendency for the x 2 to drop again toward zero 
gluino mass although in this case the x 2 /F)oF does not fall below 1 outside of the region 
of Eq. (|T4]). The x 2 /F)oF < 2 region is defined by gluino masses less than 1.2 GeV. The fit 
to the six vector quarkonia assuming the central values of Eq. (|TJ) is shown in Fig. 6. All 
three values of a 3 (M z ) picked out by the quarkonia data with light gluinos are in excellent 
agreement with the world average value for this quantity: 



a 3 (M z ) = 0.113 ± 0.003 (World Average) 0. 
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(15) 



There is at present no well established theory of supersymmetry breaking that would 
allow the unambiguous prediction of the gluino mass. Nevertheless, in the most realistic 
models that have been extensively studied, supersymmetry is softly broken, triggered by a 
super Higgs mechanism in the hidden sector of some minimal N = 1 supergravity theories 
H . Therefore, the possible SUSY breaking scenarios in such models can be parameterized 
in terms of only a few constants at the unification scale: the common gaugino mass m^2, 
scalar mass mo, and the A and B parameters of dimension mass characterizing the cubic 
and quadratic soft-breaking terms that often exist as well. Although the analysis we 
presented above is independent of any specific SUSY breaking models, it is tempting to 
discuss the implication of our results to such models. For simplicity, we now consider such 
a SUSY GUTs model which assumes the low energy form of the minimal supersymmetric 
extension of the standard model (MSSM) ||. For our purpose, it is enough to consider 
three soft-breaking parameters m^, too and A. It is interesting to note that, in this 
framework, the low gluino masses favored by our analysis are natural if the dominant 
SUSY breaking seed is the universal scalar mass mo, i.e., mi/ 2 <C too- Such a SUSY 
breaking pattern has been supported by recent considerations of proton stability in the 



context of (non-flipped) SU(5) supergravity fllPp , and also favored by cosmological studies 
TTL In fact, our results suggest 



TO1/2 = 0. (16) 

Such a model might be theoretically appealing since then SUSY breaking, like electroweak 
breaking, finds its origin in the scalar sector. On the other hand, SUSY breaking scenarios 
with mi/ 2 7^ would lead to quite large gluino masses |T!| . 

In the mi/2 = scenario, generally, one expects a supersymmetric spectrum relatively 
light compared to what one would get in other scenarios. Besides the soft- breaking pa- 
rameters and the currently unknown top quark mass, one needs two more parameters in 
order to specify the full spectrum: the ratio of the two Higgs vev's tan/3 = v%/vi, and the 
Higgs mixing parameter \i. In fact, all the gaugino masses then vanish at the tree level and 
the gauginos only receive masses through radiative corrections which are naturally small 
though dependent upon the masses of other particles In particular, the one-loop 



corrections to the gluino and photino (which is an exact mass eigenstate in this scenario) 
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masses are known, with the dominant contribution coming from graphs in which the top 
quark and its two superpartners circulate around the loop 
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where £ 2 are the masses of the two scalar top quarks (m^ < toj 2 ), and 



y) = sin26> 



.7: 



lnx — 



y 



Iny 



(17) 
(18) 

(19) 



.1 — x 1 — y 

The actual gluino and photino masses up to one-loop order are given by the absolute values 



of these mass corrections. The 9 in (^9|) is the mixing angle between the scalar partners of 
left- and right-handed top quarks, which rotates the tL,R states into the mass eigenstates 
ti 5 2- The overall factor sin2# was omitted in Ref. JT3|], corresponding to the case where 
the difference between two diagonal terms of the 2x2 mass-squared matrix of the scalar 
top quarks can be neglected (see Eq. 33 below). The more general result presented here 
has been recently calculated by one of us ||14| , and makes transparent the fact that these 
one-loop mass corrections vanish exactly if there is no left-right mixing, even if and tp, 
(now mass eigenstates themselves) have non-degenerate masses. Since we are primarily 
interested here in the gluino mass, we neglect in Eq. ([TJ) a contribution to the photino 
mass from the VF^-chargino loop diagrams [^3[. A discussion of this contribution in its 
general form will be given elsewhere [13]. As shown below, with currently favored values 
for the top quark mass mt and the two stop quark masses , m^ 2 , the function F is such 
that gluino masses below 1 GeV are quite natural. 

We now discuss the allowed regions for the relevant parameters in the mi/2 = sce- 
nario. To simplify our approach we neglect the Yukawa couplings for all the fermions 
although, strictly speaking, this approximation is only very good for the first two gener- 
ations. As a result, the diagonal elements of the sfermion mass-squared matrix can be 



written as 15 



■in. 



f 



ml + m 2 f + Mf cos2/3 



■3 J 



e/sm 2 6w 



(20) 



Here we have included the D-term contributions as well. The average mass-squared of the 
sfermions is seen to differ from the average mass-squared of the fermions by the parameter 



rn 



0' 



m 2 ^j = ml + (mj) . 



[21] 



The average mass- squared of the SUSY partners is approximately equal to the effective 
SUSY scale M| that enters into grand unification considerations. In Ref. |3J (Eq. (3.11a)) 
it was shown that the assumption of minimal SUSY unification with a light gluino (below 
Mz/2) requires the approximate relation 



M S = 150 GeV X e -518.5(sin 2 e w -0.2336) e 1.85(a 3 ^ 1 (M z )-0.113- 1 ) 



(22) 



Equating Mj with y 71 ^ an d substituting the world average values of sin 9w and olj,{Mz) 
from Ref. yields within errors the following range for mo, 



75 GeV < m < 270 GeV 



(23) 



For the top quark mass, we will assume [22 



92 GeV < m t < 147 GeV 



(24) 



In the mi/2 = scenario, the tree-level masses of the two charginos, Xi 2? are given by 

1 



m 2 ± 
x?,2 2 L 



2M Z V + =f + 4M^ 2 + 4M^cos 2 2/3 



(25) 



From Eq. (|25| ) and the requirement that the lightest chargino (xi ) nas to be heavier than 
about half the Z° mass it is found that tan/3 is restricted from both sides, i.e. 



0.441 < tan/3 < 2.266. 



(26) 



And for each value of tan/3 in the above range there is an upper limit on the absolute value 
of the Higgs mixing parameter \i. Furthermore, the lower limit of 41 GeV on the mass of 



the light CP-even Higgs boson yields the additional constraints |DJ 



0.55 < tan/3 < 0.65; or 1.5 < tan/3. 



(27) 



Combining (|26| ) with ( |27D yields two allowed ranges for the cos2/3 factor of Eq. fl20|) , 



-0.674 < cos2/3 < -0.385; or 0.406 < cos2/3 < 0.536. 



(28) 



The bounds of Eq. ( |27| ) change with the experimental lower limit on the Higgs mass, 
becoming inconsistent with Eq. (^6|) if this mass is required to be above 70 GeV. The 
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allowed parameter space of the m 1 / 2 = model also requires that the mass of the lightest 
chargino be below Mw- 

The mass-squares of the two scalar top quarks entering into Eqs. (|I?D and ( [TED are 
given by 



2 ir 

mr r = 

tl ,12 



m 2 LL + m 2 RR T \j {m\ L - m 2 RR ) 2 + Am\ R \ (29) 
with the diagonal elements (see Eq. (|2T)D) 



1 2 



m LL = mo + m 2 + M|cos2/3 - - -sm 2 6 w , (30) 



m RR = tjiq + m 2 + M§cos2/3 -sin 2 6w , (31) 



2 3 
2 - 2, 



3 

and off diagonal element 

m LR = m t (A t + ^hj) = rn t m Q At. (32) 

Here we have introduced the dimensionless mixing parameter A t as a useful combination 
of tan/3, \i and the low energy top soft- breaking parameter A t . If there is a single source 
of SUSY breaking corresponding to a single scale mo, we might expect values of A t to be 
either zero or of the order unity. In terms of At the mixing angle factor in Eq. (|T9"D is then 



sin20 = -2momA (33) 

V (™ 2 ll - m2 RR ) 2 + 4m 2 m 2 A 2 

Imposing the experimental constraint > Mz/2 together with (|23|), ( p4|) and ( |2^ ) 
requires that A t < 3.2. We also use for a 3 the value at the top mass as(m t ) — 0.1 in 
evaluating the gluino mass according to Eq. (|TTD. In Fig. 7 we show the range of gluino 
masses predicted by Eq. fll?]) as a function of A t for the allowed range of values of At, 
mo, mt, and cos2/3 described above. The constraint from Eq. (|9|) that the stop quarks 
be above half the Z° mass is also required in this allowed range. For each value of At 
there exist an upper and lower limit on the gluino mass m§. Values of At near zero are 
consistent with the near zero gluino masses of Eq. ( |136| ). Values of A t near unity are 
consistent with the x 2 minima of Eqs. (|13q| ) and (0). The entire range of values of A t 
assuming Eq. ( |23| ) yields gluino masses below 1.3 GeV in agreement with the quarkonia 
data at the two-standard-deviation level. Assuming the result of Eqs. ( |T3| ) or ([bl]) , detailed 
predictions for the individual squark and slepton masses within narrow ranges can be made. 
Much of the allowed parameter space in Fig. 7 predicts one or more of the scalar quarks 
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and leptons to have a mass between Mz/2 and Mz- If the effect of a light gluino and 
a possible light squark is taken into account the anomalously large quoted values of as 
coming from the Z° hadronic decay can perhaps be reconciled with the world average value 
Iflgfl . In addition, using a(m t ) ~ 1/127.9 in Eq. (|T8| ) the photino mass would then be 
about five times smaller than the gluino mass rrig. A stable photino of mass about 100 eV 
could provide enough dark matter to close the universe |TI| . More massive photinos would 



"overdose" the universe unless they could decay efficiently into photon plus gravitino or 
annihilate efficiently into photons. Such ultralight photinos have also been discussed as 
the explanation of other astrophysical observations ||20|1 . On the other hand photinos of 



such mass have been found to be disfavored [|21] by data from Supernova 1987A unless the 
squark masses are outside the preferred range of 60 GeV to 2.5 TeV. 

CONCLUSIONS 

We have shown that the quarkonia data behave as if there were a gluino octet in the 
region below 1 GeV. Treating the data in either of two ways, such a light gluino is favored 
by at least several standard deviations over the best QCD fits without a light gluino. In 
addition the fits without a light gluino are in conflict with SUSY unification and in dis- 
agreement with the world average values of a^iM-z) as discussed in Ref. JTJ. However, one 
must remain aware of the usual possibility that any phenomenological fit to data could be 
coincidental. We can not rule out the possibility that the real explanation for the behavior 
of the quarkonia data might lie in relativistic binding corrections or other non-perturbative 
effects, although this would contradict the general assumption that the narrowness of the 
vector quarkonia is due to perturbative QCD and asymptotic freedom. Confirmation from 
other independent data will certainly be required before the results presented here could 
be considered compelling. We are presently investigating the possibility that supporting 
evidence may be present elsewhere in existing data or that definitive experimental tests 
can be proposed. 

In addition we should address the question as to whether such light gluinos are ruled 
out by current bounds. The strongest bounds on the masses of supersymmetric particles 
come from the decay of the Z°. Any particle with electroweak charge must lie above about 
one half the Z° mass. 

Bounds on other particles such as gluinos or bounds from other processes are all to 
a greater extent model dependent. The status of these bounds is discussed in Ref. 0. 
Although many lower bounds on the gluino mass above the mass region indicated here 
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have been quoted, all of these are to some extent dependent on untested assumptions. For 
example, the stringent bounds from the hadron colliders on heavy gluinos allow windows 



for light gluinos below 50GeV. The low energy windows are well illustrated in Ref. | 23[j . 
Ref. || confirms the lack of unanimous opinion concerning whether or not very light 
gluinos have been ruled out. Many of the purported bounds have obvious loopholes some 
of which are pointed out in Ref. [|I] and elsewhere $M . For example the CUSB P5[] bound 
that disfavors gluino masses between 0.6 GeV and 2.2 GeV from the non-observation of 
7+gluinoball final states in T decay is strongly dependent on the value of the wave function 



at the origin of the gluinoball for which only models can be made [p6[ . Although the gluino 
behaves like a quark with a different color charge, it does have quartic couplings to gluons 
and gluinos that do not affect quarkonia in the same way. The binding of gluinos into new 
hadrons is therefore more closely related to the binding of gluons into new hadrons which 
is a very poorly understood area of hadronic physics at present. Cosmological constraints 
on light photinos and gluinos are subject to similar uncertainties. Most of the range 
of photino masses between 100 eV and 2 GeV is disfavored by one or more cosmological 



arguments. However there is a window noted in Ref. |?7| for a photino in the mass range 
from 4 to 15MeV. The gluino would then be in the range from 20 to 75MeV consistent 
with Eq. ( |13fo| ). It is not clear whether there is sufficient uncertainty in the cosmological 
arguments to stretch this range by a factor of four to accommodate the results of our 
Eq. fll3aD or (U). 

Given the current situation we feel that prudence requires the adoption of a conser- 
vative, non-dogmatic attitude concerning the compatibility of light gluinos with existing 
data. 
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Figure Captions 



Fig. 1. Best fit to the 3-gluon decays of six quarkonia states assuming a gluino mass of 
400 GeV. The experimental errors are increased according to "method 2" . 

Fig. 2. x 2 /6 as a function of as(Mz) for the 3-gluon decays of six quarkonia states 
treated according to method 2 but assuming a heavy gluino (mass 400 GeV). 

Fig. 3. x 2 contours for the fit to 5 vector quarkonia according to "method 1" as a function 
of as(Mz) and the gluino mass trig. 

Fig. 4. Fit to the 5 vector quarkonia using the "best fit" values from Eq. (|13a|) . 

Fig. 5. x 2 contours for the fit to 6 vector quarkonia according to "method 2" as a function 
of as(Mz) and the gluino mass m g . 

Fig. 6. Best fit to the 6 vector quarkonia with the expanded errors of method 2. The fit 
parameters are given by the central values of Eq. (14). 

Fig. 7. The allowed values of the gluino mass in the soft SUSY breaking picture with 
m i/2 = are bounded by the closed figure shown here. The range of allowed 
gluino masses at given A t corresponds to taking the full allowed range of the 
remaining parameters mo, m t , and tan/3. 
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